Preprint typeset in JHEP style - PAPER VERSION 



November 2010 



Efficient Tree-Amplitudes in Af = 4: 
Automatic BCFW Recursion in Mathematica 



Jacob L. Bourjaily 

Department of Physics, Princeton University, Princeton, NJ 08544, an d 
CN) School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540 

> 

o 

Abstract: We describe an efficient implementation of the BCFW recursion rela- 
tions for tree-amplitudes in M = 4 super Yang-Mills, which can generate analytic 
formulae for general N fc MHV colour-ordered helicity-amplitudes — which, in partic- 
ular, includes all those of non-supersymmetric Yang-Mills. This note accompanies 
the public release of the Mathematica package bef w, which can quickly (and au- 
Q_! tomatically) generate these amplitudes in a form that should be easy to export to 

any computational framework of interest, or which can be evaluated directly within 
Mathematica given external states specified by four-momenta, spinor-helicity vari- 
ables or momentum-twistors. Moreover, bef w is able to solve the BCFW recursion 

^ relations using any one of a three-parameter family of recursive 'schemes,' leading 

to an extremely wide variety of distinct analytic representations of any particular 

CN) amplitude. This flexibility is made possible by bef w's use of the momentum-twistor 

y— i Grassmannian integral to describe all tree amplitudes; and this flexibility is accompa- 

nied by a remarkable increase in efficiency, leading to formulae that can be evaluated 

t—i much faster — often by several orders of magnitude — than those previously derived 

> using BCFW. 
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1. Introduction 



The on-shell recursion relations for scattering amplitudes described by Britto, Cac- 
hazo, Feng and Witten (BCFW) [1,2] are very well known and have been widely used 
to compute scattering amplitudes for both purely-theoretical and extremely practi- 
cal purposes in a wide variety of theories [3]. They represent one of the major new 
tools in the study of quantum field theory. Theoretically, the power and simplicity of 
the recursive definitions of scattering amplitudes has allowed for the development of 
an arguably 'phenomenological' approach to the advancement of our understanding 
of quantum field theory: by making once intractable problems essentially effortless, 
many new questions can be asked — and answered. And practically, tree-amplitudes 
for processes involving many external particles are of importance for the accurate 
prediction of backgrounds for new physics at the LHC, for example; BCFW — along 
with a variety of other computational frameworks such as those based on the power- 
ful Berends-Giele recursion relations [4] — has greatly aided this effort. Considering 
for example that colour-stripped tree-amplitudes in M = 4 encode all the data of 
scattering amplitudes in ordinary, non-supersymmetric massless QCD [5], it is clear 
that understanding M = 4 is an important step along the way to understanding 
QFT in general, and as it is observed in the Standard Model as backgrounds for new 
physics at the LHC. 

Partly because of the existence and incredible simplicity of recursive definitions 
of the S-Matrix, tree-amplitudes in M = 4 have been largely understood in the liter- 
ature for some time now. Indeed, there exists today a large number of independent 
presentations of all perturbative tree-amplitudes in M = 4, including those based 
on the BCFW recursion relations [7-9], twistor string theory [10-14], contour in- 
tegrals in the Grassmannian [15,16], and the CSW recursion relations [17-19], for 
example. Many of these results were made possible in part through the existence 
of privately-developed, powerful computational tools which have proven themselves 
essential for gaining intuition and necessary for checking results. Recently, some of 
these tools have become publicly available through the release of the Mathematica 
package Gluon-Gluino-Trees (GGT), [5], which is capable of analytically computing 
all N fe MHV tree-amplitudes involving combinations of external gluons and gluinos, 
and can compute these numerically using the package £ S@M,' [6]. 

With this note, we extend the reach of these resources to include all N fc MHV 
tree-amplitudes — including those involving squarks — by making available the Math- 
ematica package bcfw, included with the submission of this paper on the arXiv. 
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In addition to its complete generality, there are two principle features of bcf w that 
should make it particularly useful to those who are interested in gaining intuition 
about or evaluating tree- amplitudes in Af = 4. First, the analytic formulae gen- 
erated by bcfw are often dramatically more compact and easier to evaluate than 
any existing formulae obtained using BCFW. To highlight the magnitude of this 
improvement, Table 1 lists the times required to evaluate unpolarized n-gluon scat- 
tering cross-sections using bcfw and GGT/SSM. 1 As will be discussed at greater length 
in section 5, these gains in efficiency can be traced directly to bcfw's: 1. use of 
momentum-twistor variables, and 2. representation of all tree-amplitudes in a fully- 
supersymmetric way (realized as contour integrals over the Grassmannian) , making 
any n-point N fc MHV helicity-amplitude easily obtained from any other. Another 
feature of bcfw that should make it useful to researchers is its ability to solve the 
BCFW-recursions using a wide variety of different recursive 'schemes,' leading to a 
large number of independent analytic formulae for any particular amplitude. 2 And 
it may be worth mentioning that the bcfw package has been designed with hopes 
of being intuitive-enough to be useful even to those with very little experience with 
Mathematica. 





mean time (10 3 s) per 
helicity component 


total time (10 3 s) 




bcfw 


GGT/SOM 


bcfw 


GGT/SOM 


A {gg ->■ ggg) 


0.11 


0.43 


2.2 


8.7 


A {gg ->■ gggg) 


0.12 


7.5 


6.1 


370 


A {gg ->■ ggggg) 


0.14 


30 


16 


3,300 


A {gg ->■ gggggg) 


0.21 


970 


49 


230,000 


A {gg -»■ ggggggg) 


0.39 


7,300 


190 


3,600,000 


A 10 (gg ->■ ggg ggggg) 


1.1 


>1,300,000 


1,100 


>1,300,000,000 


Au(gg ->• ggggggggg) 


3.1 


? 


6,700 


? 



Table 1: Evaluation-times for un-polarized n-gluon scattering cross-sections. Where in- 
dicated, estimated times are based on extrapolation from particular helicity-amplitudes. 



1 The times quoted in Table 1 (and elsewhere in this note) were obtained using the author's 
Apple laptop computer, which has a 2.6 GHz Intel Core 2 Duo processor. 

2 For example, we have included as a worked example in the demonstration hie included with 
the bcfw package the construction of all 74 linearly-independent, 20-term formulae for the 8-point 
N 2 MHV tree-amplitude, involving a total of 176 different Yangian-invariant objects. 
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Figure 1: A, 



(3) 



, +, +, +). The split-helicity 6-point NMHV amplitude. 



In[l] : = 
In [2] : = 

Out[l] : = 
Out [2] : = 



Amp [m , m , m , p , p , p] //ni ce 

Amp [m,m,m,p,p,p] //toSpinorHelicity [6] //nice 

(12) 3 (23) 3 (3451) 3 



+ 



(12) 3 (23) 3 (3561) 



(34)(4 5)(5 6)(61)(12 3 4)(2 34 5)(4 512)(512 3) (3 4) (4 5) (5 6) (6 2 3 5)(2 3 5 6)(5 6 1 2)(6 1 2 3) 

(23) 2 (34)(l|x 6 3X 34 |5) 3 (12)(2 3)(3|x 25 x 56 |l) 3 

(4 5) 3 (5 6)(61)(5|x 4 iXi2|3)s 23 S34S234 (34)(45)(61) 2 (5|x4iXi 2 |3) S6 iSi 2S612 



Figure 2: A^ ^— , — , ip^yl, +, +, ^+1/2) • ^ 6-point NMHV amplitude involving two 
gluinos and four gluons. 



in[i] := Amp[m,m,iii/2,p,p,p/2] //nice 

(12) 3 (23) 2 (3451) 2 



Qut[l] : = 



(12) 3 (23) 2 (3561) : 



(3 4)(45)(5 6)(12 34)(512 3)(4512) (3 4) (4 5) (5 6)(1 2 3 5) (5 6 1 2) (6 1 2 3) 



One of the functions defined by bcf w is 'Amp,' which can generate analytic formu- 
lae for any helicity-amplitude in M = 4. An example of how Amp can be used is given 
in Figure l. 3 Using 'm' and 'p' to denote each minus-helicity and plus-helicity gluon, 
respectively, Amp will generate any purely gluonic N fc MHV amplitude. For amplitudes 
involving 2 gluinos together with any number of gluons, a similar, simplified notation 
can be used, 4 where 'm/2' and 'p/2' indicate the two gluinos; an example of this is 
given in Figure 2. (The reader will notice that — unless 'toSpinorHelicity [n] ' is 
used — the only two kinematical invariants used by bcf w are the momentum-twistor 
'four-bracket' (••••) and its associated 'two-bracket' ( • ■ ); these will be reviewed 
along with the spinor-helicity invariants in section 2.) 



field 


SU4 .R-charge 


short-notation 


9+ 

# } 

9- 


{} 

{i} 

{i,j} 

{i,j,k} 

{1,2,3,4} 


P 

P /2(<^ {4}) 

m/2(<{=> {1,2,3}) 
m 



Table 2: Conventions for the arguments of the functions Amp, nAmp, nAmpTerms, etc. 



3 Also used in these examples is the function 'nice' which formats formulae generated by bcf w 
to be more readable — for example, by converting 'ab [1 , 2] 'h>'(1 2)'. 

4 An overall sign for these amplitudes has been implicitly fixed by the convention that the particle 
labelled 'm/2' has SU 4 i?-charge (123); refer to Table 2. 
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An example 8-point N 2 MHV helicity-amplitude involving 6 gluinos and 2 squarks. 



In[l] : = 
Out[l] : = 


Amp[{l},{l},{l},{l,3},{2,3,4},{2,3,4},{2,3,4},{2 4}] ; 
%//twistorSimplif y//nice 

(56) 2 (6 7) 2 (1236)(2345) 


(81)(12 6 7)(2 3 5 6)(2 3 6 7)(34 5 6) 


Fieure 4- ^ (5) (ib W ib {l) ^ (123) ^ (234) </> (234) </> (234) ^ {234) ib {i) ^1 

r igure 4. 1 tp + y 2 i ^+1/2 > ^+1/2 » ^+1/2' ^-1/2' ^-1/2 > ^-1/2' ^-l/2> V-l/2' ^+1/2 J • 

An example 10-point N 3 MHV helicity-amplitude involving only gluinos. 


In[l] : = 
Out[l] : = 


Amp[{l} ) {l},{l},{l},{l ) 2 ) 3},{2,3,4},{2,3,4},{2 ) 3,4} ) {2,3,4},{4}] ; 
%//twistorSimplif y//nice 

(5 6) (6 7) 2 (7 8) 2 (8 9) 2 (1 2 3 9) (2 3 4 8) (3 4 5 7) 


(10 1)(12 8 9)(2 3 78)(2 3 8 9)(3 46 7)(3 4 78)(45 6 7) 



For amplitudes involving more than two gluinos (or any number of squarks), 
simple labels such as 'm' or 'p/2' are not sufficiently precise. This is remedied by 
choosing instead to label each external particle by its SU^ .R-charge, where each of 
the external superfields are decomposed according to 



$ + = g+ + ViW+1/2 + ViVj<P {lj) + ViVjVk^l{% +ViV2V3V4 9- ■ (i-i) 

The syntactical rules which follow from these conventions are summarized in Table 2, 
but we hope they are sufficiently intuitive to be clear by example. Examples of 
how these more general helicity-component amplitudes can be specified are given 
in Figure 3, which shows an 8-point N 2 MHV helicity-amplitude involving 6 gluinos 
and 2 squarks, and Figure 4, which shows a 10-point N 3 MHV amplitude involving 10 
gluinos. These examples also illustrate the general-purpose function 'twistorSimplif y,' 
which can often greatly simplify momentum-twistor formulae. 

This paper is outlined as follows. In the next section we will review the kine- 
matics of momentum-twistors and their connection to ordinary four-momenta and 
spinor-helicity variables. In section 3, we review the tree-level BCFW recursion- 
relations as a statement about contour integrals in the momentum-twistor Grass- 
mannian, [20,21], and describe a three-parameter family of recursive 'schemes' in 
which the BCFW recursion relations can be implemented. In section 4 we describe 
the basic use of the bcf w package along with its principle functions. (A more thor- 
ough walk-through, containing numerous example computations, can be found in 
the MATHEMATICA notebook bcfw-vO-walk-through.nb distributed alongside the 
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bcf w package — attached to the submission file for this note to the arXiv.) In section 
5 we briefly discuss bcf w in the context of other existing computational tools includ- 
ing the MATHEMATICA package Gluon-Gluino- Trees (GGT), [5]. In appendix A we 
include an index of the key functions which are made available by the package bcf w. 

2. Kinematics: Momenta to Momentum-Twistors (and Back) 

By default, all tree-amplitudes generated by the bcf w package are handled internally 
as purely-holomorphic functions of the momentum-twistor variables {Z a } introduced 
by Andrew Hodges in [22] , together with an overall MHV-amplitude pre-factor which 
also depends on what is known as the 'infinity (bi-)twistor,' 1^, which associates 
with each momentum-twistor Z a a Lorentz spinor Af =1 ' 2 in the fundamental rep- 
resentation of 5X 2 (C). In addition to the many theoretical advantages of working 
with momentum-twistors, there are many indications that tree amplitudes are most 
compactly-written and most efficiently-evaluated in terms of momentum-twistors. 
But before we review this relatively novel formalism, we should reiterate that bcf w is 
fully-equipped to work with kinematics specified in terms of four-momenta or spinor- 
helicity variables (or momentum-twistors, of course), and can convert momentum- 
twistor formulae into those involving spinor-helicity variables and dual coordinates 
(but at a substantial cost in efficiency). Because of this, bcfw should be relatively 
easy to incorporate into other computational frameworks. 

The connection between ordinary four-momenta p M and momentum-twistors starts 
with the association of a (Hermitian) matrix p-- with each (real) four-momentum p^ 1 , 

p^ P ^=p»af=( Po+P3Pl - ip2 ). (2.1) 

** \Pl+ip 2 Po ~P3 J 

Noticing that p^p^ = det(p--), it follows that light-like momenta are represented 
by matrices with vanishing determinant. Any such matrix can be written as an 
outer-product, 

det(p^) = ^ = A-A-, (2.2) 

where A and A are the famous spinor-helicity variables. For real momenta, it is easy 
to see that A- = ± (A-)*, where the sign is determined by whether p M has positive 
or negative energy, respectively. Of course, this identification is only defined up-to 
an arbitrary phase: A >->■ e l6 X, A >->■ e~ %e \. Such re-phasing is induced by the action 
of little-group for massless particles in four-dimensions. 
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One of the principle advantages to working with spinor-helicity variables is that 
any function built out of the SX 2 (C)-invariants 

(XaXb) 

and [A A 6 ] 

will automatically be Lorentz-invariant up to little-group re-phasing. Amplitudes 
involving massless particles, therefore, when written in terms of spinor-helicity vari- 
ables, will be functions with uniform weight under A a h-> u\ a (with weight equal to 
minus twice the helicity of particle a). 

The next step along the road from momenta to momentum-twistors are dual 
coordinates x~ (also known as region momenta) defined (implicitly) through the 
identification 

Pa = x a - X a -i. (2.4) 

(Whenever it is necessary to fix a convention, we will choose X\ to be the origin 
of dual coordinate space.) One of the most important recent discoveries regarding 
scattering amplitudes in M = 4 SYM is that, after diving by the n-point MHV tree- 
amplitude, scattering amplitudes in M = 4 are not just superconformally-invariant 
in ordinary spacetime, but are also superconformally-invariant with respect to these 
dual-coordinates, [23,24], and this is made manifest term-by-term in BCFW, [7]. 
The existence of a conformal symmetry on this dual space led Andrew Hodges to 
propose in [22] that amplitudes be described in the twistor-space associated with 
these dual coordinates; the twistor space of dual-coordinates is known as momentum 
twistor space. 

Figure 5: The map connecting momentum-twistor variables and dual-coordinates. 




(ab) = det(A a A fe ) 
[ab] = det(A a A&) 



A a A b 



a! a 



(2.3) 



a ' 'b 
2 

a ' '6 
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Because each pair of consecutive dual coordinates are null-separated (the mo- 
menta being on-shell), the null-line joining them corresponds to a single momentum- 
twistor. And because the pair of dual coordinates (x a , x -i) encode the null-momentum 
p a , it is natural to call the momentum-twistor associated with this null-line L Z a \ 
Making this identification will associate the line (Z a ,Z a -i) in momentum-twistor 
space with the point x _i, and the line (Z a+1 , Z a ) with the point x a ; that these two 
lines intersect at the twistor Z a reflects the fact that the points x a and x a _i are 
null-separated. 

Using the conventions just established, we canonically associate a momentum- 
twistor Z a to each momentum p a according to the rule, 



Pa A a A a %a %a—l 



Z a = 



a\a 



(2.5) 



Notice that our convention of choosing x\ as the origin of dual-coordinate space triv- 



ially fixes Moreover, because this implies that p 2 = A 2 A 2 = x 2 — x 1 = x 2 , 

we see that fi 2 = £ 2c f Af(oc (A 2 A 2 }) = (° J as wei1 - Working out the rest of this map 
explicitly — as was described in [18] — we find that we may write 



^a = {Q F 1 ) ab h, where (Q/) 



ab 








/o 



(2 3) 
(2 4) (3 4) 













(25) (35) (45) 



o\ 









(2.6) 



\0 (2n) (3n) (An) • • • (n - 1 n) j 

(Q^ 1 )^ is so-named because it is a 'Formal-inverse' of the (singular) map Q ab which 
relates the /i's to the A's according to A a = Q ab jJ>b where 



Qab = 



(2n) 
(nl>(12> 

<T2> 





<nl> 



(12) 
(31) 
(12)(2 3> 




m 







<nl> 





(4 2> 









(237 


(2 3>(3 4> 











(5 3> 









(3 4>(4 5> 


(T5> '■ 



\ 



1 



n-ln) 
ln-1) 



(2.7) 



(n— In) (n-ln)(nl) 
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It is worth emphasizing that although Q a6 is singular, our conventions ensure that 
Ma = {Qf 1 ) abQbcVc and A a = Q a6 (Q^ 1 ) 6c A c , which justifies calling (Q^ 1 )^ the 
'inverse' of Q a6 . 

What we have described so far have been ordinary (Bosonic) momentum twistors; 
these have a natural extension to momentum-s'upertwistors defined by 



(2.8) 



where the Fermionic 77-components of the supertwistors are related to the ordinary 
Fermionic parameters rj which define each superfield (1.1) in precisely the same way 
that the fi variables are related to the A variables. To summarize, the components 
of the momentum-supertwistors are related to the ordinary spinor-helicity variables 
via 




AT 1 ' 2 = ZJ' 2 , and 
Aa = QabMfe, and 
Va = QabVb, and 



d=l,2 = Z 3,4 



Ma 



(Qf 1 ) 



ab 



rjb- 



(2.9) 
(2.10) 
(2.11) 



Just as spinor-helicity variables went a long way toward trivializing Lorentz- 
invariance, momentum-twistors essentially trivialize momentum conservation and 
dual conformal invariance. Momentum conservation is trivial because any set of n 
(ordered) momentum twistors will define n null-separated region momenta through 
the maps given above. Furthermore, up to little-group rescaling, dual- conformal 
transformations act on momentum-twistors as SX 4 (C) transformations, meaning 
that any function of the (only) natural SX4(C)-invariant product — namely, 'det'— 
will automatically be dual-conformally invariant if it has appropriate little-group 
weights. This suggests the natural generalization of the 'angle-bracket' (ab) defined 
for 2-spinors above would be the momentum-twistor four-bracket (••••) defined 
according to 



ab[a, b, c, d] 



(abed) = 



zl zl zl z\ 
zl zl zl z\ 
zl zl zl zl 
zt z£ z\ z\ 



Det[Zs[[{a,b,c,d}]]]; (2.12) 



So it would appear that, including also the MHV- amplitude pre-factor, all am- 
plitudes can be written in terms of four-brackets (••••) and two-brackets (• •); 
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but it is easy to see that the latter is just a special-case of the former. Notice that 
the map connecting a momentum-twistor Z a and ordinary spinor-helicity variables, 
equation (2.5), is a component-wise definition. Because any such definition is man- 
ifestly not SX 4 (C)-invariant, this map breaks dual-conformal invariance. We can 
make this clear by choosing to write Too explicitly, defining two-brackets via, 



Because momentum twistors are still somewhat unfamiliar to many researchers, 
we should mention that there is a completely canonical map between four-brackets 
and ordinary spinor-helicity variables which follows directly from definition (2.5). 
Rather than giving this map for a completely general four-bracket, we will see in the 
next section that tree-level BCFW only generates formulae involving four-brackets 
which involve at least one pair of adjacent momentum- twistors — that is, tree ampli- 
tudes involve only four-brackets of the form (a j j+1 b). Using (2.5), it is easy to 
see that 



where we have used the notation x a b = x b — x a . 5 This further simplifies in the special 
case of a four-bracket involving two pairs of adjacent momentum- twistors, 



It is worth mentioning that the fact that tree-level BCFW involves only four- 
brackets of the form (a j j+1 b) means that in general, every superamplitude in 




(a j j+1 b) = (j+1 j)(a\x aj x jb \b), 



(2.14) 



(a-l a b 6+1) = (a-1 a) (b b+1) (p a + p a+1 + ...+ p b _i + p b ) 2 
= (a-l a) (b b+l)s a ... b = (a-l a) (b 6+l)x^_ 1 b . 



(2.15) 



J\f = 4 involves strictly fewer than (M kinematical invariants 



5 This notation (and sign-convention) becomes clearer if x a b is viewed as the vector from x a to x b . 



3. Tree- Amplitudes as Contour Integrals in the Grassmannian 

The bcf w package describes each n-point N fc MHV tree-amplitude as a contour inte- 
gral in the Grassmannian G(k,n) of Axplanes in n-dimensions (see [9,15,20,25]), 



, { m= k+ 2) = 1 / d^D aa Y[ k a=l 5^{D aa Z a ) 
vo\{GL k ) J (l---k)(2---k+l)---(n---k- 



fc)(2..-A;+l).-.(n---A;-l)' 

rn,m 

= V ) 1 / d^ k D aa Yt =1 5^{D aa Z a ) ■ 

^ i vo\(GL k ) J (l---jfe)(2---fc+l) •••(«••• fc-1) i ' 

I -Da a - (dMatr ix 7 ) | =e 

= I ( residue 7) Jj5 0|4 ^(dMatrix 7 ) aa ?7 a j I , 

7er n>m L a=i J 

where we have used the scripted ',0^ m ' ) ' to indicate that this is the tree-amplitude 
divided by the (super symmetric) n-point MHV- amplitude 6 , 

n (12)(2 3)---(n-ln)(nl)' 1 ; 

As all the terms generated by the BCFW recursion relations are Yangian-invariant 
[26], they are each residues of the integral (3.1), [27,28] — computed for contours which 
'encircle' isolated poles in the Grassmannian . Therefore, each term can be described 
as a part of the complete 'tree-contour' r n m . This helps to explain the nomenclature 
of bcfw, where each superamplitude stored as a function called 'treeContour.' No- 
tice that the coefficients appearing in the Fermionic (^-functions of (3.1), dMatrix 7 , 
directly represent the isolated points in G(k,n) where the integral (3.1) develops 
a pole (of the appropriate co-dimension) which is to be 'encircled' by the contour 
r„ m , each giving rise to a particular residue of the integral. Of course, knowing the 
poles — that is, knowing just the list of points in G(k, n) (and the orientation of the 
contour about each) — is sufficient to calculate each residue using the contour integral 
(3.1); but it turns out that this is in fact unnecessary for our purposes: the BCFW 
recursion relations directly calculate the residues themselves in a canonical way. 



6 Hcre, we arc not including the ordinary momentum-conserving 5-function, <5 4 (A a A a ), because 
all momentum- twistor amplitudes are automatically on its support. 
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As described in [9], when expressed in terms of momentum-twistor variables, the 
tree-level BCFW recursion relations become the following. 



m=k+2) 



*n-l 



(3.3) 



n L ,m L 
n R ,m R 



where, 7 



Zj+i = (J+1 j) D (n-l n 1) = Z j+1 + Z 3 



n-1 n 1) 



and 



R[abcde] = 



n-lnlj) 

% = urn n (n-m i) ^ z, + z w ^ , 

Z n = nn-1 f = Z n + £ ra _i — — 

(1 J J+1 



S°^(ri a (bcde) + rj^cdea) + r] c (deab) + r]d(eabc) + r] e (abcd) 



(ab cd) (b cd e) (c d e a) (d e ab) (e ab c) 
This tree-level BCFW-bridge is illustrated in Figure 6. 

j 7 + 1 



1' X n 



(3.4) 



(3.5) 




Figure 6: The momentum-twistor BCFW-bridge (without any rotations). 

The shifted momentum-twistors in (3.4) should be understood supersymmetri- 
cally, and the shifted Fermionic ^-variables result in a shifted matrix of coefficients. 
Specifically, for terms bridged in the recursion, the residues (evaluated with shifted 
arguments) are simply multiplied, and the supersymmetric 5°' 4 's combine according 



7 It is worth noting that Zj+i and Zj are projectively equivalent; the reason for distinguishing 
them as in (3.4) is to preserve canonical little-group assignments. 
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to: 



; L ; 



a l,n L ' 



( d*i 



BCFW 
v 
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(77-1 77 lj) ' ^ (77-1 77 1 (1 jj+ln-1)' 

Thus, the tree-level BCFW recursion relations amount to little more than cutting- 
and-pasting (and re-labeling) matrices, allowing most amplitudes of interest to be 
recursed in essentially real-time. 



3.1 Generalized BCFW Recursion Schemes 



Although the recursive BCFW formula (3.3) fixes given all amplitudes with 

strictly fewer particles, (3.3) by itself does not uniquely identify any particular sum 
of residues. The reason for this is simple (and completely trivial): the lower-point 
amplitudes appearing in the recursion (3.3) can be written in any way whatsoever — 
with many choices corresponding to all the representatives r n m of each tree-contours' 
homology- class. Said another way, in order to use (3.3) to obtain a particular con- 
tour for the n-point amplitude, it is necessary to know the particular, representative 
contours for all lower-point amplitudes; but these lower-point contours need-not have 
been recursed in any particular way. In order to obtain an explicit, representative 
contour through the use of the BCFW recursion relations — i.e. using (3.3) — it is 
necessary to give a prescription for how all lower-point amplitudes are also to be 
recursed. 

One especially natural prescription would be to recurse all lower-point am- 
plitudes exactly according to equation (3.3) — with each 77-point amplitude having 



and CNvi^Hv- (3-6) 
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ordered-arguments (1, . . . , n). This is the default recursive scheme used by bcf w and 
is obtained with the function treeContour [n,m] =generalTreeContour [0,0,0] [n,m]. 
This scheme follows from Figure 6 where each lower-point amplitude is recursed pre- 
cisely according to Figure 6. 

Among the many recursive prescriptions one could imagine, a remarkable de- 
gree of complexity results from simply allowing for arbitrary (and separate) 'rota- 
tions' of the amplitudes appearing on the left- and right-hand sides of the BCFW 
bridge, 8 and also allowing for an over-all rotation of the the n-point amplitude be- 
ing recursed — or equivalently, which legs are deformed in the recursion. Specif- 
ically, letting g denote a cyclic-rotation of (an explicit formula) of an amplitude 
g : sf n {l, . . . , n) i— > £/ n (2, . . . , n, 1); then the class of generalized BCFW recursion 
schemes implemented in bcf w is given by, 

generalTreeContour[a,b,c][n,m] :g- c [^ m) ] = 9^1-1] + E ^KM ® ^K^L 

,C > n L ,m L BCFW 

n R ,m R 

where, as with the default contour prescription, this same recursive rule is used for 
every lower-point amplitude. This is illustrated in Figure 7. 




generalTreeContour[a, b, c] 



Figure 7: An illustration of the generalized BCFW recursion-schemes used by bcfw's 
function generalTreeContour [a,b, c] . Here, the legs being deformed in the left-hand 
amplitude, for example, should be thought-of as being actively 'rotated' clockwise by an 
amount 'a' relative to the default recursive scheme. 



8 When making these rotations, the homogeneous term in the recursion, ^ ra _i, must be consid- 
ered an amplitude occurring on the left. 
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By varying the parameters {a,b,c}, one can obtain a very wide- array of ana- 
lytic formulae for any particular helicity amplitude. It could be that more general 
recursion-schemes will eventually prove useful, 9 but as far as we have been able to 
check, this class of recursion schemes has proven in some sense exhaustive. Specif- 
ically, we have checked that for up to 9-particles, this three-parameter family is 
sufficient to generate all linearly independent representations of superamplitudes. 
For example, there turn out to be 74 linearly-independent formulae for the 8-point 
N 2 MHV tree amplitude, involving 176 Yangian-invariants. All of these formulae are 
worked-out explicitly as part of the demonstration file for the bcf w package. 

There are three principle reasons why researchers may find this broad-class of 
tree-amplitude formulae useful. First, knowing the range of possible tree-amplitude 
formulae helps one build intuition about amplitudes in general, and allows one to 
separate general properties about amplitudes from the peculiarities of particular 
formulae. Secondly, having many different representations available frees one from 
using unnecessarily inefficient representations of particular helicity amplitudes. For 
example, it is sometimes heard that "the" BCFW-formula (with the default scheme 
implicit) for the split-helicity amplitude is maximally-concise 10 (meaning that a max- 
imal number of terms in the tree-contour vanish); however, fixing a recursive scheme, 
this is true for at most one particular split-helicity amplitude — the other split-helicity 
amplitudes including some for which almost none of the BCFW terms vanish. And 
so, it should be possible to use the variety of representations that can be generated 
by bcf w to find a 'best-case' formula for any particular helicity amplitude of interest. 
And finally, because the BCFW formulae obtained using different recursive schemes 
often have very few spurious poles in common, it may be possible to combine a 
variety of BCFW formulae to avoid encountering spurious poles while generating 
Monte-Carlo events for phase-space integration, for example. 

It may be helpful to know that the particular recursive-scheme used by Drum- 
mond and Henn to solve the BCFW recursion relations in [7], corresponds to 
generalTreeContour [-1 , -1 , -1] ; this scheme is illustrated in Figure 8. 

9 For example, one could consider recursive schemes which make use of the parity-conjugate 
version of the BCFW-bridge, which make use of reflected (as well as rotated) lower-point amplitudes, 
or which allow for rotations of lower-point amplitudes to vary as a function of recursive depth. None 
of these generalizations are necessary for n < 9, and we suspect that this is true generally. 

1 This observation is true for the default recursion-scheme used by bcf w; in particular, the helicity 
component s^n"^ (—,...,—, + ..., +) of generalTreeContour [0,0,0] [n,m] is the gluonic ampli- 
tude with the fewest number of non-vanishing BCFW terms; but this feature is observed for very 
few of the more general recursive schemes. 
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Figure 8: Examples of particular recursion schemes, highlighting how the lower-point 
amplitudes are rotated. 



bcf w's default scheme Drummond & Henri's scheme 




generalTreeContourfO, 0, 0] generalTreeContour[— 1, — 1, — 1] 

= treeContour 



3.2 Extracting Helicity Component- Amplitudes from Tree-Contours 

To compute a particular helicity amplitude from the supersymmetric contour in- 
tegral, one need only project-out the desired Grassmann components, as dictated 
by the definition of the external superfields given in equation (1.1). Of course, 
the component fields of are given in terms of r/-variables, which, as described in 
section 2, are related to the momentum-supertwistor Grassmann parameters r] a via 

Va=(Q- F 1 ) ab r ]b - (3-7) 

Because the matrix of coefficients of the Grassmann 77's for each residue is nothing 
but its corresponding dMatrix, we have that 

D aa rj a = D ab (Qp^baVa = C aa rj a . (3.8) 

In terms of the Grassmannian integral (3.1), this means that we may write 

k k 

(residue) <5 ' 4 ((dMatrix) aa ?7 a ) = (residue) <5 ' 4 ((dMatrix) Q , 6 (QabInverse[n]) 6a r/ a ) 

a=l a=l 

k 

= (residue) JJ 5° |4 ((cHatMatrix) aa r/ a ) . (3.9) 
0=1 
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Upon explicitly including the full MHV super-amplitude we obtain, 



fresidue 



(12>..>1> 
fresidue 



I! 6 ° 14 i^aVa) I] 5 °' 4 ((cHatMatrix)^) 



a=l 
fc+2 



a=l 



(3.10) 



(12>...(nl> 



}J<5 0|4 (cMatrix Sa ^), 



3=1 



where we have defined the matrix Caa according to 





( Cu 


Cl2 ' 


' Cln-1 




r*~ — 


Cfci 


Cfc2 ' 


" Ckn-1 


Cfc n 






4 ■ 


' A"_! 


As 




U 


Al ■ 


A n-1 


X 2 




(3.11) 



It is worth noting that just as each dMatrix represents an isolated point in the 
Grassmannian of A;-planes in n-dimensions, each cMatrix gives an isolated point in 
the Grassmannian of m{= k + 2)-planes in n-dimensions. Indeed, these are the iso- 
lated poles 'encircled' by the (original) twistor-space Grassmannian contour-integral 
of [15], 



,( m =fc+2) 



vol(GL 



?L m ) j (1 



m)(2- • -ra+1) 



In- 



m-1) 



(3.12) 



The momentum-twistor Grassmannian integral (3.1) was derived from the original 
twistor-space integral (3.12) in [20], where it was shown how the MHV-prefactor 
arises naturally as the Jacobian of the change- of -variables in going from the (space- 
time) twistor-space variables W a to the momentum-twistor-space variables Z a . 

Now, having the matrix of coefficients of the rj- variables, it is particularly simple 
to extract any helicity component amplitude. For example, pure-glue amplitudes are 
given by 



•^n V ■ ■ i J\ i ■ ■ ■ i Jmi ■ ■ ■ I 



residue^) 



(3.13) 



4- (i 



7er„ 



(12). ..(nl> 



Det[cMatrix 7 [[All, . . . , j m }]]]) 
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More generally, each helicity amplitude can be 'projected-out' of the superam- 
plitude by multiplying each residue in the tree-contour by the appropriate set of four 
(m x m)-minors of its corresponding matrix Csa- The list of minors which project-out 
a particular helicity component-amplitude is given by the function parselnput [] . 

4. The bcfw MATHEMATICA Package 

A separate Mathematica notebook — distributed along with bcfw.m — has been pre- 
pared to introduce the reader to the many functions of bcfw and their primary usage. 
We hope that the demonstration notebook is sufficiently self-contained for most users. 
In this section, we briefly describe the basic algorithmic structures which underly the 
bcfw package, with an emphasis on the features that are likely to prove useful beyond 
the limited framework of Mathematica. 

4.1 Setup and Initialization 



Initialization of the package is simple: so long as the file being used has been saved to 
the same directory as the package's source bcfw.m, one need only call the following: 



In[l] : = 
Out[l] : = 


SetDirectory [NotebookDirectory [] ] ; 
<<bcfw.m 








j j+l 










7 BCFW 










1 n 








Efficient Tree-Amplitudes in =4 SYM 

via BCFW in the Momentum-Twistor Grassmannian 
Jacob L. Bourjaily, 2010 





4.2 Getting Started with Analytic Tree Amplitudes 

To start gaining intuition for how helicity-amplitudes can be specified in bcfw, con- 
sider a very simple example: the 8-point MHV amplitude A 8 (+, +, — , +, +, — , +, +). 
This amplitude can easily be found using bcfw through the command, 
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in[i]:= Amp[p,p,m,p,p,m,p,p] 

ab[3, 6] 4 



Out[l] : = 



ab[l, 2]ab[2, 3]ab[3, 4]ab[4, 5]ab[5, 6]ab[6, 7]ab[7, 8]ab[8, 1 



To make the result more aesthetically appealing, any output of bcfw can be 
wrapped by the function 'nice [] ' which formats the result so that it is more "human- 
readable." For example, using nice, the above command would return: 



in[i] := Amp [p,p,m,p,p,m,p,p] //nice 

(36> 4 



Out[l] : = 



;i 2) (2 3) (3 4) (4 5) (5 6) (6 7) (7 8) (81) 



We have chosen to make 'nice' formatting an ' opt-in' option for users, so that the 
underlying structure is transparent at all times — and in order to avoid the pitfalls 
of conditional formatting in Mathematica while maximizing the ease of symbolic 
manipulation. 

Although the analytic formulae for tree amplitudes quickly become too long and 
complex for visual comprehension, bcfw's function Amp will in fact write-out any 
amplitude. As one further example, consider the 6-point NMHV alternating helicity 
amplitude. 



in[i]:= Amp [m,p,m,p,m,p] //nice 

(15) 4 ((35)(1234) - (34)(1235)) z 



Out[l] : = 



(12)(2 3)(3 4) (4 5)(5 6)(61)(12 3 4)(12 3 5)(124 5)(134 5)(2 3 45) 
((13)(56)(1235) - (15)((3 6)(1235) + (3 5) (2 3 6 l))) 4 

(12) (2 3) (3 4) (4 5) (5 6) (6 1)(1 2 3 5)(1 2 5 6)(1 3 5 6)(2 3 5 6) (2 3 6 1) 
((13)(56)(1345) - (15)((36)(1345) + (34)(1356) + (3 5)(3 4 6 l))) 4 

(12) (2 3) (3 4) (4 5) (5 6)(6 1)(1 345)(1 35 6)(345 6)(346 1)(45 6 1) 



We should emphasize, however, that direct evaluation of the formulae generated 
by Amp (or AmpTerms) are often dramatically less efficient than what can be obtained 
using nAmp (or nAmpTerms). 11 This directly reflects the efficiency gained by the 
momentum-twistor Grassmanniannian representation of superamplitudes. 12 



11 This is true even with fairly intelligent caching. Because of this, researchers interested in 
transferring the formulae generated by bcfw to other frameworks should seriously consider using 
the superamplitudes directly. 

12 To better understand this, observe that each cMatrix includes as its first fc-rows the matrix 
cHatMatrix=dMatrix. Qablnverse [n] ; this introduces many new kincmatical invariants into each 
term — the two-brackets — while simultaneously duplicating each column of dMatrix many times, 
greatly obfuscating an underlying simplicity with fundamentally redundant information. 
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Table 3: 6-point NMHV superamplitude Ag> , given by treeContour [6,3]. 





Name 


residue 


dMatrix 


1. 


#[12345] 


1 


((2345) (3451) (4512) (5123) (1234) 0) 


(12 34)(2 3 45)(34 51)(4512)(512 3}(12 34) 


2. 


#[13456] 


1 


((3456) (4561) (5613) (6134) (1345)) 


(13 45)(3 4 5 6)(45 61)(5 613)(613 4)(13 45) 


3. 


#[1235 6] 


1 


((2356) (3561) (5612) (6123) (1235)) 


(1 2 3 5) (2 3 5 6) (3 5 6 1) (5 6 1 2) (6 1 2 3} (1 2 3 5} 



As described in the previous section, each superamplitude is represented by 
bcfw as a contour integral in the momentum-twistor Grassmannian (3.1). The 
particular representation of the n-particle N^^MHV superamplitude derived via 
the BCFW recursion scheme with rotations {a,b,c} is obtained with the function 
generalTreeContour [a,b,c] [n,m] (see section 3.1). The default representation — 
obtained using the default recursion scheme, with {a,b, c}={0 , 0,0}, is obtained with 
treeContour [n,m] . For example, the default representation of the 6-point NMHV 
superamplitude is given in Table 3. 

4.3 Referencing, Generating, or Specifying Kinematical Data 

In order to evaluate amplitudes numerically using bcfw, kinematical data must first 
be defined. This can be done by calling upon a list of reference momentum-twistors, 
freshly-generating random kinematics, or by specifying kinematical data explicitly: 

1. useRef erences [n] : use a standard set of reference momentum-twistors; these 
reference twistors were carefully selected so that 

• all components are integer- valued (and small); 

• there are no physical or spurious singularities; 

• all kinematical invariants are uniformly positive (that is, s a .„b > for all 
ranges a ... b), and that these invariants are numerically given by ratios of 
relatively small integers — leading to amplitudes that are ratios of integers 
that are 'not-too-horrendously-long'; 
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Table 4: Reference momentum-twistors used in bcfw's function useRef erences [n] . 
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Figure 9: Evaluation of 10-point N 3 MHV helicity amplitudes to infinite precision using 
reference momentum-twistors. The timing reflects the fact that the first computation 
determined the full superamplitude and projected-out a particular helicity component, 
while the second only needed to perform the projection. 



useRef erences [10] ; 

nAmp [m,m,m,m,m,p,p,p,p,p]/ /withTiming 

Evaluation of the 10-point N 3 MHV amplitude required 46.7. ms to complete. 

17886892256634020134576330754470391777 
280278666971743564282064966167680000 
nAmp [m,p,m,p,m,p,m,p,m,p]/ /withTiming 

Evaluation of the 10-point N 3 MHV amplitude required 8.6. ms to complete. 

5007045380847632725336670465304701314367799201604575059832902148541 
213450466354689126392301641566350924968168379805192061706240000 



In Table 4 we give a sample of the reference momentum-twistors. Because 
an arbitrary set of momentum-twistors define on-shell, momentum- conserving 
kinematics, there are no constraints from momentum conservation. Therefore, 
choosing simply the first n twistors from the list in Table 4 will suffice. It 
is worth mentioning, however, that these reference momentum-twistors are 
neither canonically normalized 13 , nor do they map to real four-momenta in IR 3,1 . 

Nonetheless, reference twistors are extremely well-suited for debugging, 
checking identities, and finding relations to infinite precision. As one can see 
in Figure 9, using bcfw's built-in reference momentum-twistors can quickly 
lead to scattering amplitudes that are known to infinite-precision. Notice 
that in Figure 9, once the superamplitude has been computed for any helicity- 
component, all subsequent components are obtained quite rapidly. 

13 By not having canonical normalization, we mean that there are non-trivial, Lorentz-frame (and 
hence also little-group)-dependent kinematical scale-factors in the spinors; however, this tends to 
only cause a problem when combining/comparing multiple helicity component-amplitudes. 



In[l] : = 

Dut[l] : = 
In [2] : = 



Figure 10: Evaluation of 12-point N 4 MHV helicity amplitudes with random kinematics. 





useRandomKinematics [12] ; 


In[l] : = 


nAmp [m,p,m,p,m,p,m,p,m,p,m,p]/ /withTiming 




Evaluation of the 12-point N 4 MHV amplitude required 596. ms to complete. 


Qut[l] : = 


-274.127 - 5171.81 I 


In [2] : = 


nAmp [p,m,p,m,p,m,p,m,p,m,p,m]/ /withTiming 




Evaluation of the 12-point N 4 MHV amplitude required 71.4. ms to complete. 


Out [2] : = 


-274.127 + 5171.81 I 



2. useRandomKinematics [n] : use randomly-generated kinematics in IR 3 ' 1 . This 
function chooses a random set of (optionally rational or arbitrary-precision) on- 
shell four-momenta in IR 3 ' 1 , and sets up essentially all the kinematical variables 
of potential interest, including 

• momentum-twistors {~^} =Zs, given as an (rax 4) matrix — the n rows list- 
ing the four homogeneous components of each momentum-twistor; 
useRandomKinematics [n] also defines the 'dual' momentum-twistors 
{W} =Ws, which, although not used by bcfw, may be found useful by 
some researchers; 

-> ^ 

• spinor-helicity variables { A } =Ls and { A } =Lbs, each an (n x 2) ma- 
trix of components; these have been normalized so that A a = ± (A a )*, as 
described in section 2; 

• fourMomenta, an (n x 4) matrix of the components (p° , p x , p y , p z ) of each 
four- momentum; 

• regionMomenta, the dual-coordinates (described in section 2), given as a 
n-length list of 2 x 2 Hermitian matrices; 

An example of using random kinematics is shown in Figure 10, where the 
two alternating-helicity 12-point N 4 MHV amplitudes were evaluated. Notice 
as before that once the superamplitude has been evaluated, subsequent helicity 
components are quickly extracted. Also, observe that the randomly-generated 
spinors and momentum-twistors have been appropriately normalized so that 
parity-conjugation results in complex-conjugation of the amplitude. 
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3. using user-defined kinematics, given in terms of: 

(a) setupUsingFourMomenta[f ourMomentaList] : generates momentum-twistors 
and spinor-helicity variables for the input list of four-momenta, 

f ourMomentaList, which must be given as an n-tuple of four- vectors list- 
ing the components of each four-momentum; the list of four-momenta 
must conserve momentum; 

(b) setupUsingSpinors [Ls ,Lbs] : generates momentum-twistors given the 

— > ~^ 

spinor-helicity variables Ls= { A } and Lbs= { A } each given as an (n x 2) 
matrix of components; 

(c) setupUsingTwistors [twistorList] : establishes the necessary kinemat- 
ical functions given the (unconstrained) list of user-generated momentum- 
twistors. 

Examples of how each of these functions can be used can be found in the 
demonstration file included with the bcf w package. 

4.4 Numerical Evaluation of Tree Amplitudes 

As has been emphasized throughout this paper, the principle sources of bcf w's ef- 
ficiency are manifest supersymmetry and the use of momentum-twistor variables, 
which are both made manifest in the momentum-twistor Grassmannian integral (3.1). 
Because these ingredients — or at least their implementation — are quite novel in bcf w, 
it is worth describing in some detail how amplitudes are evaluated numerically by 
the bcf w package. 

The basic evaluation strategy is outlined in Table 5, where we give the basic 
evaluation times for each step in the evaluation of the 10-point N 3 MHV alternating- 
helicity tree-amplitude. 

Because of the central role played by momentum-twistors, the first step of any 
numerical evaluation is the establishment of momentum-twistor variables which can 
then be used to compute the kinematical invariants that determine any scattering am- 
plitude. This can be done in a number of different ways — as described in the previous 
subsection. Although this should be completely clear from the discussions above, this 
step is not very computationally-intensive (and indeed, can be discounted entirely 
by choosing to randomly-generate momentum-twistors instead of four-momenta). 

Because of the ubiquity of the MHV- amplitude pre-factor, 1/((1 2) • • • (n 1)), and 
the map (Q^ 1 ) , used to relate the momentum-twistors' ^-variables to the 77 variables 
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Table 5: The general evaluation strategy used by bcf w, with a break-down of evaluation- 
time requirements for each step in the case of the alternating-helicity 10-point N 3 MHV 
tree- amplitude, A^j (— , +, — , +, — , +, — , +, — , +) (for random kinematics). 

1. setupUsingRandomKinematics [10] 

(a) generate random (on-shell, rational, momentum-conserving) 

four- momenta in R 3 ' 1 ; define spinors and momentum-twistors 3.61 ms 

(b) evaluate the universal objects nMHVpref actor and nQinverse 1.34 ms 

2 . nAmp [m,p,m,p,ni,p,m,p,ni,p] 

(a) evaluate the full-superamplitude, which is stored stored as the 



of the external superfields, bcf w evaluates these two objects and stores them globally 
whenever new kinematical data is defined. 

The first step in the evaluation of any particular helicity amplitude is actually 
the evaluation of the full superamplitude — represented as the list of BCFW-terms, 
where each is described by the pair {residue , dMatrix} (which is stored in memory 
as the function nContour [a,b, c] [n,m]). Because particular helicity amplitudes 
are usually specified with respect to the 77-variables of the external superfields, the 
dMatrix of each residue is then converted to the corresponding cMatrix as described 
in section 2. 

Once each BCFW-term has been evaluated numerically and stored as the pair 
{residue , cMatrix}, it is relatively easy to extract any particular helicity component 
amplitude — by multiplying each residue by the appropriate four (m x m) minors 
of its corresponding cMatrix. This last step is nothing exotic: it is merely the 
evaluation of the Grassmann integrals J YliLi d°\ 4 rji which project-out a helicity- 
component amplitude from the superamplitude. 




23.2 ms 
3.03 ms 
4.01 ms 



Total Time: 



35.2 ms 
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4.5 Example Applications 

In the demonstration file which accompanies the bcf w package, several examples are 
given which illustrate how bcf w can be used as a tool to verify results, find identities, 
or learn about scattering amplitudes more generally. In particular, these examples 
emphasize how using integer-valued reference momentum-twistors to compute am- 
plitudes (and individual BCFW-terms) to infinite-precision can prove quite useful 
theoretically. The examples include: 

• a verification of supersymmetric Ward identities; in particular, we check one of 
the 'cyclic' identities described in [29] for the 10-point N 3 MHV amplitude — 

n_ a(5) (Am ^(3) .(1) .(4) ,(3) ,(24) ,(14) ,(12) ,(234) (1234) \ 
u_ "^10 I <P -1/2' r+1/2' y+ 1/2' Y+ 1/2' r + 1/2' V^O ' 9o 5 Y>0 '^-1/2' »- J 

' ^lO I ^-1/2' ^+1/2' ^+1/2' ^+1/2' ^+1/2' V^O 'V^O ) ¥>0 iW-\/2i9- J 

+ A {5) (*b im *b (4) ^ (4) ^ (3) ^ {23) <^ (14) <^ (12) ^ (234) o {1234) I 

-|-^ 10 ^-1/2' ™+l/2' ™+l/2' ™+l/2' ™+l/2' ™0 'V^O ' V^O 'r_l/2'i/- J 

+ ^ (5) ^ (123) ^ (4) ^ (3) ^ (3) ^ {24) ^ (14) ^ (12) ^ (234) o {1234) ^l 

"1-^10 I ^i^) f+1/2) f+1/2' r + l/2i % ' *P0 ' ^0 5</ ; _l/2'y- J 

—the verification of which is illustrated in Figure 11, highlighting the power of 
knowing amplitudes to infinite precision; 

• an explicit verification of the t/i-decoupling identity for the 10-point N 3 MHV 
tree-amplitude (which, although a trivial consequence of any Lagrangian field 
theory, is a highly non-trivial check of numerical code! 14 ); 

• a complete classification of the linearly-independent BCFW-generated formulae 
for the 8-point N 2 MHV supersymmetric tree-amplitude. 



We thank Freddy Cachazo for this suggestion. 
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Figure 11: Verifying a supersymmetric Ward identity of the 10-point N 3 MHV amplitude. 



useRef erences [10] 

m[i] := List [nAmp[{l, 2, 3}, {3}, {1}, {4}, {3}, {2, 4}, {1, 4}, {1, 2}, {2, 3, 4}, {1,2, 3, 4}] 

nAmp[{l, 2, 3}, {4}, {1}, {4}, {3}, {2, 4}, {1, 3}, {1, 2}, {2, 3, 4}, {1,2, 3, 4}] 
nAmp[{l, 2, 3}, {4}, {1}, {4}, {3}, {2, 3}, {1, 4}, {1, 2}, {2, 3, 4}, {1,2, 3, 4}] 

nAmp[{l, 2, 3}, {4}, {1}, {3}, {3}, {2, 4}, {1, 4}, {1, 2}, {2, 3, 4}, {1,2, 3, 4}] 

79370862801471295255 
28753113503920775424' 

1275513453387873135869428633786428491 



0ut[l] : = 



77923676342112832490222204964602880 
40428898488502522106856665437052838463 

~ 10951273590541549612279689882333035520 ' 

16319258699414773847825256760953737) 

~ 1057119835135513498965174929610240 J 

In[2]:= Total[0ut[l]] 
Out [2] : = 



5. Comparison of bcf w with Other Computational Tools 

Although we do hope that the bcf w package will prove useful in computations directly 
relevant to collider physics, its primary role will likely be as a tool for gaining intuition 
about scattering amplitudes, checking results/conjectures, and as a working example 
of a novel computational strategy that could perhaps be implemented much more 
efficiently by researchers with more computational expertise by optimizing either 
hardware or software. 

Although there exists a wide variety of tools for computing scattering amplitudes — 
including COMIX [30], AMEGIC++ [31], CompHEP [32], MadGraph [33], HELAC [34], and 
ALPHA [35] — it would be difficult for us to make any just comparison between these 
packages and bcf w. (But it would be very interesting to see how the representation 
of amplitudes used by bcf w would compare with the results of [36], or the impressive 
algorithms described in [37] based on the Berends-Giele recursion relations [4].) This 
is both because of the inherent inefficiencies of any M athematica package relative 
to compiled code, and also because these packages make use of a wide variety of 
different computational strategies and have widely-different scopes of purpose. 

Nevertheless, at least within the limited scope of M athematica, the recent 
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release of the public package Gluon-Gluiono-Trees (GGT) [5] has made at least 
a passing comparison between the two packages justified. After all, both pack- 
ages were written in and for Mathematica, both are based on the BCFW re- 
cursion relations, and both are (in principle) capable of computing all pure-glue 
amplitudes in Af = 4, giving ample room for comparison. Moreover, because GGT 
uses the form of tree-amplitudes in Af = 4 obtained by Drummond and Henn 
in [7] — which corresponds to the recursion scheme with {a,b, c}={-1 , -1 , -1} (see 
section 3.1), the amplitudes computed in GGT match term-by-term the output of 
nAmpTerms [-1,-1,-1] [helicityConf iguration] , it is natural to wonder how the 
two packages compare in efficiency (at least in gross terms). 

The two packages GGT and bcfw were compared in passing in section 1, where 
we listed in Table 1, the times required to compute unpolarized (colour-ordered) n- 
gluon scattering cross sections. Considering that both packages were written for 
largely-theoretical purposes (as demonstrations of algorithms, for intuition-building, 
and checking results), it is easy to argue that the computation of all 1002 helicity- 
amplitudes for 10 gluons is not the computation either package was designed for. 

A more reasonable comparison, and one which may shed light on the source of the 
disparity between the two packages, would be the evaluation of an individual helicity 
amplitudes. Consider for example the 8-particle N 2 MHV alternating-helicity ampli- 
tude. This computation is illustrated in Figure 12. Here, we have explicitly separated 
the time required by GGT to convert the analytic amplitude to a form suitable for nu- 
merical evaluation and the time actually required for evaluation by S@M. 15 Of course, 
only the latter time is of essential interest in computations, as this represents (perhaps 
poorly) the essential complexity of the analytic forms of the amplitudes generated by 
the two packages' frameworks. (We should mention in passing that Figure 12 is not 
really a fair comparison between the two packages: in its 4 milliseconds, bcfw actu- 
ally evaluated the entire 8-particle N 2 MHV superamplitude — only a half-millisecond 
of which was used to project-out the alternating-helicity component-amplitude.) 

A better understanding the relative efficiency between the two packages requires 
a more systematic survey than the isolated example Figure 12. To give an idea of how 
the two packages compare more generally, Table 6 lists the kinematically-averaged 
evaluation times for a range of particular pure-glue scattering amplitudes. Similar 

15 As indicated in Figure 12, over half of the computation time is absorbed by the third-party 
package used by GGT, 'S@M' [6], to evaluate formulae written in its spinor-helicity formalism. Al- 
though this package is far from optimal, it seems unlikely to account for more than a small fraction 
of the disparity between the two packages. 
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Figure 12: Comparing the evaluation of A 8 (+, — , +, — , +, — , +, — ) in GGT/SOM vs. bcfw. 



In[l] : = 


spmorr orm=GGI toSpmors LCiGInnmnvgluon |_o , 1 ,4 , oj J //witnTiming; 


In [2] : = 


N [spinorForm] / /withTiming 




Rva 1 nat i nn o"F p f nnrt i on iT-C-r'TtnSTH n nr s ypmii ypri A OQ cjpr nnti r 1~ n rrmrnl pt" p 




Evaluation of the function N required 5.84 seconds to complete. 


Out [2] : = 


-0.395021 + 0.3107191 


In [3] : = 


nAmp [p , m , p , m , p , m , p , m] / /withTiming 




Evaluation of the 8-point N 2 MHV amplitude required 4.30 ms to complete. 




"4 4) , 9- ,9+,9-,9+,9-,9+,9-) ■ 


Out [3] : = 


-0.395021 + 0.3107191 



time comparisons could have been made for amplitudes involving gluinos (it is worth 
remembering that because bcfw always computes the entire superamplitude before 
extracting a component, the times quoted in Table 6 would be roughly the same 
regardless of the helicity components used for comparison). 

The scale of the differences observed in Table 6 is hard to overlook, and naturally 
raises the question of what underlies the difference in efficiency? This question 
seems especially relevant considering that if the particular BCFW recursion scheme 
nAmpTerms [-1 , -1 , -1] is used to compute amplitudes in bcfw the two packages agree 
term-by-term for every component amplitude (including those involving gluinos). We 
suspect that the two main sources of relative efficiency between bcfw and GGT are 
simply: the use of momentum-twistor variables, and keeping supersymmetry manifest 
throughout every computation, by describing amplitudes directly as contour integrals 
in the Grassmannian. 

A good illustration of the relative simplicity afforded by momentum-twistor vari- 
ables is the comparison between bcfw and GGT when using Drummond and Henn's 
recursive scheme, so that both packages are computing essentially the same functions. 
Indeed, merely translating the formulae for amplitudes given in [7] into momentum 
twistors would seem to offer a remarkable improvement. This has in fact been done 
for NMHV [15, 21] and N 2 MHV [38] tree-amplitudes, but not more generally (al- 
though we suspect it would not be very difficult). 

To see how merely translating the spinor-helicity formulae written by Drum- 
mond and Henn could result in a remarkable increase inefficiency, observe that when 
written in terms of spinor-helicity variables and dual coordinates, each term in an 
N fc MHV amplitude will typically involve (see, e.g. [7]) a large number of general- 
ized spinor-helicity angle-brackets of the form ' (ai\x aia2 x a2a3 ■ ■ ■ x akak+1 \a,k+i)\ with 
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Table 6: Comparing evaluation-times in bcfw and GGT/S@M for particular helicity amplitudes. 





Time (1(T 3 s) 








NMHV Amnlitndes 




GGT 




Time (10~ 3 s) 


a( 3 )( | | ^ 


0.67 


16 


N 2 MHV Amnlitndps; 


bcfw 


GGT 


1 1 u) 


0.98 


71 




0.91 


230 


4 3) (-+,- +,-+,+) 


1.4 


240 


4 4) (- +,-,+, -,+,-) 


2.0 


2,300 


4 3) (- +,-,+, -,+,+,+) 


1.9 


690 


4 4) (- +,- +, -,+,-,+) 


4.2 


16,000 


4 3) (-+>- +>-+>••• ,+) 


2.6 


1,100 


4 4) (- +,-,+, -,+,-,+,+) 


8.6 


72,000 


4 3 o(- +,-,+,-,+,••• ,+) 


3.5 


2,000 


4 4 o(-, +,-,+,-,+,-, +,+,+) 


16 


260,000 


4?(- +,-,+, -,+,••• ,+) 


4.6 


3,400 


4 4 1 ) (- +,-,+, -+,-+,... ,+) 


30 


740,000 


4 3) (- +,-,+,-,+,••• ,+) 


5.7 


5,000 




48 


1,900,000 





Time (10~ 3 s) 




N 3 MHV Amplitudes 


bcfw 


GGT 




A? ( , +, , +, , , ) 


1.1 


2,500 




"^■8 ( 5 "I" 5 5 + 5 ? + ? 5 ) 


3.5 


97,000 




"4-9 ( _ J "I") _ ) "I") _ ) ~K _ J "1"? _ ) 


12 


1,100,000 




"4-10 ( — ) ~K ~-> ~K ~) +5 ~-> +5 — ) +) 


35 


14,000,000 




"4-11 ( _ J +5 _ ) +) _ ) +) ~) +) _ ) +! +) 


97 


? 




( _ ) +) ~) +5 ~) +5 _ 5 +5 _ J +) +5 +) 


210 


? 



(A; + 1) region momenta sandwiched-between the spinors A ai and A afc+1 . This leads 
to an essentially boundless variety of new kinematical invariants that must be com- 
puted separately (the number of which grow very rapidly as k increases). In contrast, 
as described in section 2, the number of momentum-twistor four-brackets that can 
occur for tree-amplitudes is strictly bounded. To understand the magnitude of this 
handicap, consider that the expression used by GGT/S@M to compute the 10-point 
N 3 MHV alternating-helicity tree-amplitude involves more than 256-thousand distinct 
kinematical functions; 16 in contrast, the formulae obtained with the Drummond and 
Henn recursion scheme in the momentum-twistor Grassmannian involves only 130 
separate four-brackets. 17 This makes the difference between 35 milliseconds and 4 
hours given in Table 6 much easier to comprehend. 



16 Here, we quote the number of distinct spinor-helicity invariants that would be used by the 
package S@M for numerical evaluation. 

17 Of course, these invariants are not all independent: they are related in complicated ways through 
momentum conservation. The fact that momentum-conservation is trivial for momentum-twistor 
variables has a dramatic impact on the simplicity of the formulae that result. 
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6. Conclusions 



We have described a general, versatile, and efficient implementation of the tree-level 
BCFW recursion relations within the framework of Mathematica which has been 
realized by the bcf w package which is included with the submission of this posting 
on the arXiv. 18 

Having access to an efficient, reliable, flexible, and robust toolbox for computing 
scattering amplitudes in M = 4 has proven an essential resource, and a important 
source of theoretical 'data.' It is hard to overlook the exciting recent advances that 
have been made in our understanding of scattering amplitudes, and many of these 
results have relied heavily on being able to decisively rule-out or quickly confirm a 
wide-array of new ideas and proposals, leading to many new insights, and helping to 
establish what has a chance to become a fundamentally new descriptions of quantum 
field theory. 

We hope that the bcf w package proves itself useful to a wide range of researchers — 
both as a reliable and efficient black-box for computing amplitudes, and as an ed- 
ucational resource for gaining intuition about the still somewhat unfamiliar, but 
extremely powerful new tools available to describe amplitude such as the momentum- 
twistor Grassmannian that have played an important role in the recent development 
of our understanding of scattering amplitudes in Af = 4. 
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project's page on http://hepforge.org, where it will be generally maintained by the author. 
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A. Glossary of Functions Defined by bcf w 



A.l BCFW Recursion in the Grassmannian 



bcf wPartitions [n,m] : gives a list of the pairs {(n^, mi), {iir, m^)} which 
should be bridged-together to compute the tree-level n-point N^ rn ~ 2 ^MHV am- 
plitude An according to 



A { m h) 



E 

(n L ,m L ) BCFW 
(n R ,m R ) 



A( m n) 



generalBCFWbridge [a,b, c] [{n L , m L }, {n R , m R }] : computes the contributions 
to a given tree-amplitude arising from the term (^) A^ R \ where the 



BCFW 



arguments of An"l L \ A^^ and An have been 'rotated' relative to (1, ... , n*) 
by a, b and c, respectively. There are at least two ways to envision how these 
different recursive schemes are defined. As illustrated in the figure on the left 
below, one can view the left- and right-side amplitudes as being actively rotated 
by amounts a and b, respectively — where rotation means that the legs which 
are deformed at the next stage of the recursion are offset relative to the default 
positions. On the right-hand figure, we view the same legs as being deformed 
at every stage, but at each stage of the recursion the amplitudes on the left- 
and on the right- have been cyclically-relabeled according to the figure (which 
induces the 'rotations' of the first description). 

c + j c + O'+i) 




• generalTreeContour [a,b,c] [n,m] : the complete supersymmetric tree-amplitude 
An™^ — given in the representation obtained from the BCFW recursion scheme 
generalBCFWbridge [a, b, c] . The contour is given as a list of contributing 
terms of the form 

{residue, dMatrix}, 
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where each dMatrix represents a point in the Grassmannian G(k, n) — in the 
gauge appropriate to the form of residue. The full amplitude is given by the 
contour integral over the Grassmannian which 'encloses' each of these poles; the 
contour integral's contribution from each pole is its corresponding residue. A 
given helicity-component of each residue is obtained by multiplying the function 
residue by the appropriate minors of its dMatrix — which can also be viewed 
as the matrix of coefficients of the 77's appearing in the Fermionic ^-functions. 

• treeContour [n,m] : the default representation of supersymmetric tree-amplitudes. 
Specifically, treeContour[n, m] = generalTreeContour[0, 0, 0] [n, m]. 

A. 2 Component Amplitude Extraction 

• Amp [helicityConf iguration] : returns an analytic formula for the ampli- 
tude specified by helicityConf iguration, obtained using the default BCFW 
recursion scheme. Specifically, 

Amp[helicityConf iguration] = Total [AmpTerms[0, 0, 0] [helicityConf iguration]]. 

• AmpTerms [a_ : , b_ : , c_ : 0] [helicityConf iguration] : generates each BCFW 
term for the specified helicityConf iguration analytically, using the recur- 
sive scheme specified by a,b,c. We should emphasize that formulae given by 
AmpTerms are not particularly well-suited for numerical purposes: each term 
contains a great-deal of redundancy that would be systematically removed by 
the function nAmpTerms. 

• mhvPref actor [n] : the universal pre-factor for all amplitudes obtained within 
the momentum-twistor Grassmannian integral: 

mhvPref actor[n] : = — -. (A.l) 

1 J ab[l, 2] • • • ab[n, 1] v ' 

• parselnput [helicityConf iguration] : parses the input to the functions Amp, 
AmpTerms, nAmp, and nAmpTerms, which for an N m ~ 2 MHV amplitude returns a 
list of four m-tuples which specify the four (m x m)-minors of each cMatrix 
which project-out the desired helicity configuration. Consider for example the 

amplitude ^4 4) (^2/ 2 , ^+i/ 2 » ^+i/2> 0o 13) > ^?V2' ^-V2' ^o^) of Figure 3. 

That amplitude was computed via the command 

Amp [{!}, {!}, {!}, {1, 3}, {2, 3, 4}, {2, 3, 4}, {2, 3, 4}, {2 4}] ; (A.2) 
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the arguments of the function Amp were first parsed by parselnput, returning 



m[i]:= parselnput [{1} , {1} , {1} , {1 , 3} , {2,3 ,4} , {2 ,3,4} , {2 , 3,4} , {2 4}] 
Out[i]:= {{2, 1,3, 4}, {5, 6, 7, 8}, {4, 5, 6, 7}, {5, 6, 7, 8}} 



meaning that the specified helicity component-amplitude is obtained by multi- 
plying each residue in the treeContour by 

Det [cMatriXi[[All, {2, 1, 3, 4}]]]Det [cMatrix^fAll, {5, 6, 7, 8}]]] 
x Det [cMatriXi[[All, {4, 5, 6, 7}]]] Det [cMatrix ± [[All, {5, 6, 7, 8}]]] . ^ A ' 3 '' 

If the input is ill-formed, apparently in error, or otherwise confusing (to bcfw), 
parselnput will print an error message and return ' — 1.' 

A. 3 Input, Translation, and Random-Generation of Kinematical Data 

• Qab[n] : the matrix defined in equation (2.7) which provides the map from 
momentum-twistors' components fi a to the associated spinor-variables A a via 

Aa = QabfMb. (A.4) 



• Qablnverse [n] : the matrix defined in equation (2.6) — a "formal inverse" 
of Q ab — which provides a canonical map from spinor-helicity variables A a to 
momentum-twistor components \i a . Of course, because Q ab is a singular map, 
a particular representative of the one-to-many map (Q^ 1 ) has was chosen in 
order to make equation (2.5) literally correct once we have fixed the origin for 
the space of dual coordinate to be the point x\. 

• setupUsingFourMomenta[inputFourMomentaList] : translates a given list of 
momentum-conserving four-momenta (given as an n-tuple of four-vectors) into 

a standard set of spinor variables Ls= A = {A a } a=1 n and Lbs= A = {A } a= i . n 
— each, a globally-defined (n x 2) matrix — together with a set of corresponding 
momentum-twistors Zs= Z = {Z a } a =i...n- 

• setupUsingSpinors [inputLs , inputLbs] : translates a given list of momentum- 
conserving spinor variables inputLs and inputLbs (to become Ls and Lbs), 
each given as an (n x 2) matrix, into a standard set momentum-twistor variables 
Zs, given as a globally-defined (n x 4)-matrix. 
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• setupUsingTwistors [input TwistorList] : translates a given list of generic 
momentum-twistors (given as an (rax 4) matrix) into a standard set of momentum- 
conserving spinor variables Ls= A = {A a } a=1 n and Lbs= A = {A a } a= i... n 
— each, a globally-defined (ra x 2) matrix. 

• useRandomKinematics [n] : randomly-generates on-shell, momentum-conserving 
four-momenta for n-particles in Minkowski signature, and sets up essentially 
all kinematical variables of potential interest, including 

— momentum-twistors {^} =Zs, given as an (rax4)-matrix — the ra rows list- 
ing the four homogeneous components of each momentum-twistor; 
useRandomKinematics [n] also defines the 'dual' momentum-twistors 
{W} =Ws, which, although not used by bcfw, may be found useful for 
some researchers; 

_> ^ 

— spinor- helicity variables { A } =Ls and { A } =Lbs, each an (ra x 2)-matrix 
of components; these have been normalized so that A a = ± (A a )*, as de- 
scribed in section 2; 

— fourMomenta, an (ra x 4) matrix of the components (p° , p x , p y , p z ) of each 
four- momentum; 

— regionMomenta, the dual-coordinates (described in section 2), given as a 
ra-length list of 2 x 2 Hermitian matrices; 

• useRandomSpinors [n] : the same random-kinematics-engine as that behind 
useRandomKinematics, but much faster because it neglects to define any of the 
superfluous kinematical quantities (which are defined by useRandomKinematics, 
making it rather inefficient). 

• useRandomTwistors [n] : an extremely fast, minimalistic kinematics-generator 
that picks random, integer-valued twistors, and uses these to define the cor- 
responding { A } =Ls, etc. As with useRef erences [n] , the kinematical data 
generated using useRandomTwistors [n] will neither be canonically normal- 
ized, nor will they correspond to real four-momenta in Minkowski signature. 
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• useRef erences [n] : uses a standard references-set of momentum-twistors to 
define external kinematics; these reference twistors were carefully selected so 
that 

— all components are integer-valued and relatively small; 

— there are no physical or spurious singularities; 

— all kinematical invariants are uniformly positive (that is, s a ,__j, > for 
all ranges a ... b), and these invariants are numerically given by rational 
numbers composed of small integers — leading to amplitudes that are ratios 
of integers that are 'not-too-horrendously-long.' 

The first few reference momentum-twistors are listed in Table 4. Because 
an arbitrary set of momentum-twistors define on-shell, momentum-conserving 
kinematics, there are no constraints from momentum conservation. Therefore, 
choosing simply the fist n from the list in Table 4 will suffice. 

A. 4 Numerical Evaluation of Tree Amplitudes 

All of the functions below assume that the global variables such as the list Zs have 
been appropriately defined by the user as described in section A. 3. 

• nAmptTerms [a_: 0,b_: , c_: 0] [helicityConf iguration] : The numerically- 
optimized analogue of AmpTerms. This function will return the numeric expres- 
sions which contribute to the given amplitude. It has built-within it several 
optimizations which (often dramatically) improve the evaluation of a given 
super-amplitude. Researchers interested in exporting the general computa- 
tional strategy of the bcf w to more efficient computational environments should 
try to mimic the functionality of nAmpTerms rather than using the analytic for- 
mulae generated by AmpTerms. 

• nAmp [helicityConf iguration] :=Total [nAmpTerms [0,0,0] [helicityConf iguration] 

• toN: includes the replacement rules necessary to numerically evaluate each 
momentum-twistor four-bracket and associated two-bracket. If the global list 
of momentum twistors Zs have been defined, then (expression) /. toN will 
convert expression to Numbers. 
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A. 5 ^Esthetics 

• nice [expression] : formats expression by replacing ab [x- • • y] i— )■ (x---y), 
and by writing any level-zero matrices in MatrixForm. 

• niceTime [timelnSeconds] : converts a time measured in seconds timelnSeconds, 
to human-readable form. 



in[i] : = niceTime [299792458] 
niceTime [1CT 14 ] 

out[i]:= 9.51 years 
10.0 is 



• order [expression, (option)] : the the optional option option set to 1 — its 
default value — order will order all angle-brackets to a canonical ordering which 
prioritizes x-like pairing of arguments in four-brackets, and cyclically-ordered 
two-brackets (picking up any necessary signs from the necessary permutations). 
For example, 



in[i]:= order [ab [1,3, 4, 9]] 
out[i]:= -ab[3,4,9,l] 



Notice that order has guessed that the expression should be Zg-ordered (if 
there had been more brackets in expression involving particle-labels larger 
than 9, order would have chosen the maximum argument of all ab for cyclic 
ordering. 

If the optional argument option were set to 0, then order will lexicographically- 
order the arguments of all angle-brackets, picking-up all necessary signs from 
the permutations involved. 

• twistorSimplif y [expression, (option)] : a general-mess-of-a function, which 
uses order together with FullSimplif y endowed with the power of the most 
elementary, three-term Schouten identity 

(Xab)(Xcd) + (Xbc)(Xad) + (Xca)(Xbd) = (A.5) 

(where X is any bi-twistor) to simplify expressions. twistorSimplif y can 
occasionally yield favourable results — especially when expression is relatively 
simple, but — like FullSimplify — is liable to exhaust patience. 
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• withTiming [expression] : an all-purpose timing-wrapper, which will print- 
to-screen the time required to evaluate any expression using human-readable 
units of time (see niceTime). withTiming will also identify the function being 
evaluated in the output, and includes special formatting for the functions nAmp, 
Amp, nAmpTerms, etc. 
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